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Abstract. This paper presents large scale extinction 
maps of most nearby Giant Molecular Clouds of the 
Galaxy (Lupus, p Ophiuchus, Scorpius, Coalsack, Taurus, 
Chamaeleon, Musca, Corona Australis, Serpens, IC 5146, 
Vela, Orion, Monoceros Rl and R2, Rosette, Carina) de- 
rived from a star count method using an adaptive grid 
and a wavelet decomposition applied to the optical data 
provided by the USNO-Precision Measuring Machine. The 
distribution of the extinction in the clouds leads to esti- 
mate their total individual masses M and their maximum 
of extinction. I show that the relation between the mass 
contained within an iso-extinction contour and the extinc- 
tion is similar from cloud to cloud and allows the extrap- 
olation of the maximum of extinction in the range 5.7 to 
25.5 magnitudes. I found that about half of the mass is 
contained in regions where the visual extinction is smaller 
than 1 magnitude. The star count method used on large 
scale (~ 250 square degrees) is a powerful and relatively 
straightforward method to estimate the mass of molecular 
complexes. A systematic study of the all sky would lead 
to discover new clouds as I did in the Lupus complex for 
which I found a sixth cloud of about 10 4 Mq. 
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1. Introduction 

Various methods have been recently developed to estimate 
the mass of matter contained in giant molecular clouds 
(GMC) using millimetric and far infrared observations 
([Boulanger et al. 1998); [Mizuno et al. 1998Q. Neverthe- 



less, the mapping of the optical/near-infrared extinction, 
based on star counts still remain the most straightforward 
way to estimate the mass in form of dust grains. These 
maps can be usefully compared to longer wavelength emis- 
sion maps in order to derive the essential physical param- 
eters of the interstellar medium such as the gas to dust 
mass ratio, the dumpiness of the medium or the optical 
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and morphological properties of the dust grains. The star 
count method was first proposed by Wolf (1923) and has 
been applied to Schmidt plates during several decades. It 
consists to count the number of stars by interval of mag- 
nitudes (i.e. between m — 1/2 and m + 1/2) in each cell 
of a regular rectangular grid in an obscured area and to 
compare the result with the counts obtained in a suppos- 
edly unextinguis hed re gion. In order to improve the spatial 
resolution, Bok ( 1956 ) proposed to make count up to the 
completeness limiting magnitude (to < mu m ). Since the 
number of stars counted is much larger when compared 
to counts performed in an interval of 1 magnitude, the 
step of the grid can be reduced. Quite recently extinc- 
tion maps of several southern clouds have been drawn by 
Gregorio Hetem et al. ( 1988 ) using this secon d met hod. 
Even more recently, Andreazza and Vilas-Boas ( 1996 ) ob- 
tained extinction map of the Corona Australis and Lupus 
clouds. Counts were done visually using a x30 magnifi- 
cation microscope. With the digitised Schmidt plates the 
star counts method can be worked out much more eas- 
ily across much larger fields. The star count methods have 
also tremendously evolved thanks to the processing of dig- 
ital data with high capacity computers. Cambresy et al. 
(1997), for instance, have developed a counting method, 
for the DENIS data in the Chamaeleon I cloud, that takes 
advantage of this new environment. The aim of this paper 
is to apply this method to the optical plates digitised with 
the USNO-PMM in a sample of Giant Molecular Clouds 
to derive their extinction map (Vela, Fig. [I]; Carina, Fig. 
||; Musca, Fig. [| Coalsack, Fig. f|; Chamaeleon, Fig. 5; 
Corona Australis, Fig. g; IC 5146, Fig. Q Lupus, Fig. |; 
p Ophiuchus, Fig. ^|; Orion, Fig. [ll]; Taurus, Fig. [l|; Ser- 
pens, Fig. |l3|). The method is shortly described in Sect. |[ 
In Sect. ||, deduced parameters from the extinction map 
are presented and Sect. ^ deals with individual clouds. 



2. Star counts 

2.1. Method 

The star counts method is based on the comparison of lo- 
cal stellar densities. A drawback of the classical method 
is that it requests a grid step. If the step is too small, this 
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may lead to empty cells in highly extinguished regions, 
and if it is too large, it results in a low spatial resolution. 
My new approach consists in fixing the number of counted 
stars per cell rather than the step of the grid. Since un- 
certainties in star counts follow a poissonian distribution, 
they are independent of the local extinction with an adap- 
tive grid where the number of stars in each cell remains 
constant. Practically, I used a fixed number of 20 stars per 
cell and a filtering method involving a wavelet decomposi- 
tion that filters the noise. This method has been described 



in more details by Cambresy ( 1998 ) 



Fig. 1. Extinction map of Vela from B counts (J2000 co- 
ordinates) 



plates (6 < —35°) in blue and red. Internal photometry 
estimators are believed to be accurate to about 0.15 mag- 
nitude but systematic errors can reach 0.25 magnitude in 
the North and 0.5 magnitude in the South. Astrometric 
error is typically of the order of 0.25 arcsecond. This ac- 
curacy is an important parameter in order to count only 
once those stars which are detected twice because they are 
located in the overlap of two adjacent plates. 

All the extinction maps presented here have been 
drawn in greyscale with iso-extinction contours overlaid. 
On the right side of each map, a scale indicates the corre- 
spondence between colours and visual extinction, and the 
value of the contours. Stars brighter than the 4 th visual 



magnitude (Hofflcit and Jaschek 1991) are marked with a 
filled circle. 



Fig. 2. Extinction map of Carina from R counts (J2000 
coordinates) 



Fig. 3. Extinction map of Musca from R counts (J2000 
coordinates) 

I have applied the method to 24 GMCs. Counts and fil- 
tering are fully automatic. Because of the wide field (~ 250 
square degrees for Orion), the counts must be corrected for 
the variation of the background stellar density with galac- 
tic latitude. Extinction and stellar density are related by: 



Ay 



-log 

a 



Drefjb) 

D 



(1) 



where A\ is the extinction at the wavelength A, D is the 
background stellar density, D re f the density in the refer- 
ence field (depending on the galactic latitude b), and a is 
defined by : 



log(D re/ ) 



est. 



(2) 



where m\ is the magnitude at the wavelength A. 
Assuming an exponential law for the stellar density, 
D re f(b) — Z?oe~ Q ' b ', a linear correction with the galactic 
latitude b must be applied to the extinction value given 
by Eq. (|l|). The correction consists, therefore, in subtract- 
ing log[D re f(b)] = log(Do) — a\b\ log (e) to the extinction 
value A\(b). This operation corrects the slope of A\{b) 
which becomes close to zero, and set the zero point of ex- 
tinction. All maps are converted into visual magnitudes 



assuming an extinction law of Cardelli et al. (1989) for 



which At 



1.337 and 



0.751. 



Here, the USNO-PMM catalogue (Monet 1996) is used 



to derive the extinction map. It results from the digiti- 
sation of POSS (down to —35° in declination) and ESO 



2.2. Artifacts 

Bright stars can produce artifacts in the extinction maps. 
A very bright star, actually, produces a large disc on the 
plates that prevents the detection of the fainter star (for 
example, a Crux in the Coalsack or Antares in p Ophi- 
uchus, in Fig. [I] and ||, respectively). The magnitude of 
Antares is my — 0.96, and it shows up in the extinction 
map as a disc of 35' diameter which mimics an extinction 
of 8 magnitudes. Moreover Antares is accompanied of re- 
flection nebulae that prevent source extraction. In Fig. 
|ll] the well known Orion constellation is drawn over the 
map and the brighter stars appear, e Ori (Alnilam) the 
central star of the constellation, free of any reflection neb- 
ula, is represented by a disc of ~ 14' for a magnitude of 
my = 1.7. Fortunately, these artifacts can be easily iden- 
tified when the bright star is isolated. When stars are in 
the line of sight of the obscured area, the circularity of 
a small extinguished zone is just an indication, but the 
only straightforward way to rule out a doubt is to make a 
direct visual inspection of the Schmidt plate. 

Reflection nebulae are a more difficult problem to iden- 
tified since they are not always circular. Each time it was 
possible, I chose the R plate because the reflection is much 
lower in R than in B. B plates were preferred when R 
plates showed obvious important defaults (e.g. edge of the 
plate) . 



2.3. Uncertainties 

Extinction estimations suffer from intrinsic and system- 
atic errors. Intrinsic uncertainties result essentially from 
the star counts itself. The obtained distribution follows 
a Poisson law for which the parameter is precisely the 
number of stars counted in each cell, i.e. 20. Equation (|l]) 
shows that two multiplicative factors, depending on which 
colour the star counts is done, are needed to convert the 
stellar density into visual extinction. These factors are a 
(the slope of the luminosity function (|j)), and the conver- 
sion factor A\/Ay. For B and R band, the derived extinc- 
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Fig. 4. Extinction map of Coalsack from R counts (J2000 coordinates) 

Fig. 5. Extinction map of the Chamaclcon complex from B counts (J2000 coordinates) 

Fig. 6. Extinction map of Corona Australis from B counts (J2000 coordinates) 



tion accuracies are ^023 magnitudes and j^ ' 4 magnitudes, 
respectively. 

Also, for highly obscured region, densities are esti- 
mated from counts on large surfaces, typically larger than 
~ 10'. It is obvious that, in this case, the true peak of 
extinction is underestimated since we have only an av- 
erage value. The resulting effect on the extinction map 
is similar to the saturation produced by bright stars on 
Schmidt plates. This effect cannot be easily estimated and 
is highly dependant on the cloud (about ~ 80 magnitudes 



in p Ophiuchus, see Sect. 3.2) 



Moreover, systematic errors due to the determination 
of the zero point of extinction are also present. Extinction 
mappings use larger areas than the cloud itself in order to 
estimate correctly the zero point. This systematic uncer- 
tainty can be neglected in most cases. 



3. The extinction maps and the derived 
parameters 

3.1. Fractal distribution of matter in molecular clouds 

Fractals in molecular clouds characterize a geometrical 
property which is the dilatation invariance (i.e. self- 
similar fractals) of their structure. A fractal dimension in 
cloud has been first found in Earth's atmospheric clouds 
by comparing the perimeter of a cloud with the area of 
rain. Then, using Viking images, a fractal structure for 
Marti an clouds has also been found. Bazell and Desert 
obtained similar results for the interstellar cirrus 



(£98L 

discovered by IRAS. Hetem and Lepinc ( 1993 ) used this 
geometrical approach to generate clouds with some sta- 
tistical properties observed in real clouds. They showed 
that classical models of spherical clouds can be improved 
by a fractal modelisation which depends only on one or 
two free parameters. The mass spectrum of interstellar 
clouds can also be understood assuming a fractal struc- 
ture (Elmegreen and Falgarone 1996). Larson (1995) went 
further, showing that the Taurus cloud also presents a 
fractal structure in the distribution of its young stellar 
objects. 



Blitz and Williams (1997) claim, however, that clouds 
are no longer fractal since they found a characteristic size 
scale in the Taurus cloud. They showed that the Taurus 
cloud is not fractal for a size scale of 0.25-0.5 pc which may 
correspond to a transition from a turbulent outer envelope 
to an inner coherent core. This is not inconsistent with a 
fractal representation of the cloud for size scales greater 



Table 1. Cloud properties. Distances are taken from lit- 
erature, masses (expressed in solar masses) are defined by 
the regression line log M( Ay) = log(Mrot) + a x Ay, 
the maxima of extinction, Ay, is measured from star 
counts, and A v , is extrapolated from the previous equa- 
tion assuming a fractal structure and the last column is 
the value of the slope a 



Cloud Name 


d (pc) 


Am, 
A V 


4 e 

Ay 


M To t 


Slope 


Lupus I 


100 (1) 


5.3 


7.1 


10 4 


-0.56 


Lupus II 


100 (1) 


3.8 


5.7 


80 


-0.33 


Lupus III 


100 (1) 


4.9 


7.6 


1150 


-0.40 


Lupus IV 


100 (1) 


5.3 


7.0 


630 


-0.40 


Lupus V 


100 (1) 


5.2 


10.6 


2500 


-0.32 


Lupus VI 


100 


4.8 


7.0 


10 4 


-0.57 


p Ophiuchus 


120 (1) 


9.4 


25.5 


6600 


-0.15 


Scorpius 


120 


6.4 


7.0 


6000 


-0.54 


Taurus 


140W 


7.5 


15.7 


1.1 10 4 


-0.26 


Coalsack 


150 (D(3) 


6.6 


6.3 


1.4 10 4 


-0.63 


Musca 


150« 


5.7 


10.1 


550 


-0.27 


Chamaeleon III 150 (4) 


3.7 


7.8 


1300 


-0.40 


Chamaeleon I 


160 (4) 


5.2 


12.9 


1800 


-0.25 


Corona- 


170 (1) 


5.4 


10.7 


1600 


-0.30 


Australis 












Chamaeleon II 


178 (4) 


4.9 


12.3 


800 


-0.22 


Serpens 


259 (5) 


10.1 


??? 


??? 


??? 


IC 5146 


400 (6) 


6.5 


11.9 


2900 


-0.29 


Vela 


500 (7) 


4.0 


??? 


??? 


??? 


Orion 


500 (8) 


7.5 


20.3 


310 5 


-0.27 


Crossbones 


830 (8) 


4.4 


10.4 


7.3 10 4 


-0.47 


Monoceros R2 


830 (8) 


4.1 


10.5 


1.2 10 5 


-0.48 


Monoceros Rl 


1600 (9) 


5.4 


20.8 


2.7 10 5 


-0.26 


Rosette 


1600 (9) 


8.4 


20.4 


510 5 


-0.28 


Carina 


2500 (10) 


8.3 


82? 


5.5 10 5 ? 


-0.07? 


(1) : Knude 


and Hog ( 


1998), 


(2) 


: Kenyon 


at al. 



(|1994j), (3) : Fran co (099J ), (4) : Whitt et et al. ( |1997| ), 
(5) : Straiz ys et al. ( [1996Q , (6) : Lada et al (|l994| ), (7) : Dun- 
canet al. ( jl996| ), (8) : Madd alena et al. ( |l986j ), (9) : Turner 
(|1976[), (10) : Feinstein (|l995|) 



than 0.5 pc. Fractal in physics are defined over a number 
of decades and have always a lower limit . 

Fractal structures in clouds can be characterized by a 
linear relation between the radius of a circle and the mass 
that it encompasses in a log-log diagram. Several defini- 
tions of fractal exist, and this definition can be written 
M oc L , where L is the radius of the circle and D the 
fractal dimension of the cloud. The mass measured is, in 
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Fig. 7. Extinction map of IC5146 from R counts (J2000 coordinates) 

Fig. 8. Extinction map of the Lupus complex from B counts (J2000 coordinates) 

Fig. 9. Extinction map of p Ophiuchus and Scorpius clouds from R counts (J2000 coordinates) 
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Fig. 10. Mass contained inside the iso-extinction contours 
versus the extinction (solid line), and regression line for 
the linear part. Annotations indicate the area in square 
degrees contained by the iso-extinction contours Ay. 

fact, contained in a cylinder of base radius L and of unde- 
fined height H (because the cloud depth is not constant 
over the surface of the base of the cylinder). In our case, 
we are interested in the relation between the mass and the 
extinction. Since the extinction is related to the size H, 
which represents the depth of the cloud as defined above, 
we seek for a relation between mass and Ay (or H) , L be- 



ing now undefined. The logarithm of the mass is found to 



vary linearly with the extinction over a range of extinction 
magnitudes (Fig. Il0|). We have: 



logM = log M To t + slope x A\ 



(3) 



This result is compatible with a fractal structure of the 
cloud if Ay oc log if, i.e. if the density of matter follows a 
power law, which i s, precisely, what is used in modelling 
interstellar clouds ( Bernard et al. 1993 ). 



3. 2. Maximum of extinction 

Fig. |l0| shows the relation between iso-extinction contours 
and the logarithm of the mass contained inside these con- 
tours for the Taurus cloud (see extinction map in Fig. |l2| ). 
The relation is linear for Ay 5.5. For higher extinction, 
mass is deficient because the star count leads to under- 
estimating the extinction. Indeed, for highly extinguished 
regions, the low number density of stars requires a larger 
area to pick up enough stars and estimate the extinction. 
The result is therefore an average value over a large area 



in which the extinction is, in fact, greater. In the Tau- 
rus cloud, I found that this turn off occurs for a size of 



~ 1.7 pc. According to Blitz and Williams (1997) a turn 
off toward higher masses should appear for a size scale of 
~ 0.5 pc. Obviously, our value corresponds to a limitation 
of the star counts method with optical data and not to a 
real characteristic size scale of the cloud. 

It is therefore natural to extrapolate, down to a min- 
imum mass, the linear part of the relation log[M(A„)] 
versus Ay to determine a maximum of extinction. This 
maximum is obtained using the regression line (H) and 
represents the higher extinction that can be measured, 
would the cloud be fractal at all size scales. As Blitz and 



Williams ( 1997 ) have shown, there is a characteristic size 
scale above which the density profile becomes steeper. The 
extrapolation of the linear relation gives, therefore, a lower 
limit for the densest core extinction. Derived values of Ay 
are presented in the 4 th column of Table [l| and correspond 
to a minimum mass of LMq. This minimum mass is a typi- 
cal stellar mass and, an extrapolation toward lower masses 
would be meaningless. 

Except for the Carina cloud, maxima are found in the 
range from 5.7 to 25.5 magnitudes of visual extinction 
with a median value of 10.6 magnitudes. 

We stress the point that extinction can be larger. The 
p Ophiuchus cloud ,for example, is known to show extinc- 
tion peaks of about ~ 100 magnitudes (Casanova et al 



1995 ), whereas we obtain only 25.5 magnitudes. Never- 
theless, this value compared to the 9.4 magnitudes effec- 
tively measured indicates that we need deeper optical ob- 
servations or near-infrared data — such as those provided 
by the DENIS survey — to investigate more deeply the 
cloud. The Coalsack and Scorpius are the only clouds for 
which measured and extrapolated extinctions are similar 
(see Sect. |). 

3.3. Mass 

Assuming a gas to dust ratio, the mass of a cloud can be 
obtained using the relation (Dickman 1978) : 



M 



{ad) 2 fi-^-^A v {i 



A 



where a is the angular size of a pixel map, d the dis- 
tance to the cloud, p the mean molecular weight cor- 
rected for helium abundance, and i is a pixel of the extinc- 
tion map. Uncertainties on the determination of masses 
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come essentially from the distance which is always dim- 
cult to evaluate. Assuming a correct distance estimation, 
error resulting from magnitude uncertainties can be eval- 
uated: an underestimation of 0.5 magnitude of visual ex- 
tinction implies a reduction of the total mass of a fac- 
tor ~ 2. According to Savage and Mathis ( 1979| ), the gas 



to dust ratio is 4 s " = 1-87 x 10 21 cm 2 .mag 1 where 



N H = N H i + 2N H2 . Kim and Martin ( |1996| ) show that 
this value depends on the total to selective extinction ra- 
tio Ry — Av/Eb-v- For Ry = 5.3, the gas to dust ratio 
would be divided by a factor 1.2. The value of Ry is sup- 
posed to be larger in molecular clouds than in the general 
interstellar medium but the variation with the extinction 
is not clearly established. So, I used the general value of 
3.1 and the gas to dust ratio of Savage and Mathis. Mass 
of the cores of the clouds may, therefore, be overestimated 
by a factor ~ 1.2. 

In Fig. |l^, the relation log[M(A,j)] vs. Ay extrapo- 
lated toward the zero extinction gives an estimation of 
the total mass of the cloud using Eq. (||). Masses obtained 
are shown in the Table [I]. The median mass is 2900Mr, 
and th e range is from 80 to 5 10 5 M^. Using expression (|T 



side the iso-extinction curve 1.0 magnitudes. This value 
is remarkably stable from cloud to cloud with a standard 
deviation of 0.3. 

4. Remarks on individual clouds 

4-1. Vela and Serpens 

In the Vela and the Serpens clouds (Fig. [I] and [l3[ respec- 
tively), there is no linear relation between log [M (A v )) and 
Ay . Extrapolation for the maximum of extinction or for 
the total mass estimations is not possible. However, mass 
lower limits can be obtained using the extinction map di- 
rectly: 5.71O 4 M and 1.110 5 M Q for Vela and Serpens, 
respectively. It is difficult to understand why there is no 
linear relation for these two clouds, even for low values of 
extinction. 



4-2. Carina 

The study of the Carina (Fig. |2|) presents aberrant values 
for the slope of the line log [M(.A„)] (see Table [l]). Conse- 
quently, the maximum of the extrapolated extinction, 82 
magnitudes, cannot be trusted. The important reflection 
in the Carina region is probably responsible for the shape 
of the extinction map. Stars cannot be detected because of 
the reflection and thus, extinction cannot be derived from 
R star counts. Infrared data are requested to eliminate the 
contribution of the nebulae. 



4-3. Musca-Chamaeleon 

The Chamaeleon I ha s already been mappe d with DENIS 
star counts in J band (Cambresy et al. 1997) and the max- 



imum of extinction was estimated to be ~ 10 magnitudes. 
Using B star counts we find here 5.2 magnitudes. This dif- 
ference is normal since J is less sensitive to extinction than 
B. The important remark is that the extrapolated value 
for the extinction derived from B star counts is 12.9, con- 
sistent with the value obtained with J star counts. We ob- 
tain the same result for the Chamaeleon II cloud for which 
J star counts lead also to a maximum of about ^10 mag- 
nitudes whereas the extrapolated value from B counts is 
12.3 magnitudes. Infrared data are definitely necessary to 
investigate the cores of the clouds. 

Besides, the shape of the whole Musca-Chamaeleon 
extinction and the IRAS 100/im maps are very similar. 
Disregarding the far-infrared gradient produced by the 
heating by the galactic plane, there is a good match of 
the far-infrared emission and extinction contours and fila- 
mentary connections. This region is well adapted to study 
the correlation between the extinction and the far-infrared 
emission because there is no massive stars which heat the 
dust. The 100/xm flux can, therefore, be converted in a 
relatively straightforward way into a col umn density unit 



we alsc remark that half of the total mass is located out- 1998) 



using the 60/100/im colour temperature ( Boulanger et al 



4-4- Coalsack and Scorpius 

The extinction map of the Coalsack is displayed in Fig. 

The edge of the cloud contains the brightest star of 
the Southern Cross (a = 12 h 26 m 36 s , S = -63°05'57"). 
This cloud is known to be a conglomerate of dust mate- 
rial, its distance is therefore difficult to estimate. Franco 



(1995) using Stromgren photometry give s a di stance of 



150-200 pc. More recently Knude and Hog ( 1998 ) estimate 
a distance of 100-150 pc using Hipparcos data. Finally, 
I adopted an intermediate value of 150 pc to derive the 
mass of the Coalsack. Maximum extinction estimations 
for this cloud are 6.6 and 6.3 for the measured and the 
extrapolated values, respectively. The lowest limit for the 
maximum of extinction, as defined in Sect. |3.2| , is reached, 
but no characteristic size scale has been found by studying 
the shape of log [M(^4„)] at high extinction. These values 
are too close, regarding their uncertainties, to reflect any 
evidence of a clumpy structure. 

Nyman et al. ( 1989] ) have made a CO survey of the 
Coalsack cloud. They have divided the cloud into 4 re- 
gions. Regions I and II which correspond to the northern 
part of the area near a Crux, are well correlated with 
the extinction map. Nyman et al. have defined two other 
regions which have no obvious counterpart in extinction. 
Region III below -64°of declination in the western part 
of the cloud is not seen in the extinction map and this 
may result of a scanning defect of the plates. The same 
problem exists for the region IV which is a filament in the 
eastern part of the cloud at S ~ —64°. 

The Scorpius cloud is located near the p Ophiuchus 
cloud (Fig. H). The 120 pc distance used to derive its mass 
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Fig. 11. Extinction map of Orion, Monoceros I, Rosette, Monoceros II, Crossbones from R counts (J2000 coordinates) 



Fig. 12. Extinction map of Taurus from R counts (J2000 coordinates) 



Fig. 13. Extinction map of Serpens from R counts (J2000 coordinates) 



is the p Ophiuchus distance ( Knude and Hog 1998 ). As for 
the Coalsack cloud, measured and extrapolated extinction 
are similar : 6.4 and 7.0, respectively. But, no evidence of 



in these regions. The complex has been first separated into 
4 clouds ( [Schwartz 1977 ), and then, a fift h cloud has been 
recently discovered using 13 CO survey (Tachihara et al 



a characteristic size scale can be found. For both clouds, 



this result is not surprising since the maximum of mea- 
sured extinction reaches the extrapolated value. Would 
dense cores with steeper extinction profile be found, the 
measured extinction would have been significantly greater 
than the extrapolated value. 

4-5. Corona Australis 

The extinction in this cloud has recently been derived us- 
ing star counts on B plates by Andreazza and Vilas-Boas 
( |1996| ). Our methods are very similar and we obtain, there- 
fore, comparable results. The main difference comes from 
the most extinguished core. Since they use a regular grid, 
they cannot investigate cores where the mean distance 
between two stars is greater than their grid step. Con- 
sequently, they obtained a plateau where I find 4 distinct 
cores. 



4.6. IC5146 



CO observations are presented in Lada et al. ( 1994 ). The 
two eastern cores in the 13 CO map have only one coun- 
terpart in the extinction map because the bright nebula, 
Lynds 424, preven ts the star detections in that region. 
Lada et al. (1994) also present an extinction map of a 
part of the cloud derived from H — K colour excess obser- 
vations. This colour is particularly well adapted for such 
investigations because, in one hand, infrared wavelengths 
allow deeper studies and, on the other hand, H — K colour 
has a small dispersion versus the spectral type of stars. I 
obtain similar low iso-extinction contours but they reach 
a much greater maximum of visual extinction of about 20 
magnitudes. 

4-7. Lupus 



New estimations of distance using Hipparcos data (Knudc 
and Hog 199S| ) have led to locate the Lupus complex (Fig. 
S|) at only 100 pc from the Sun. Lupus turns out to be the 
most nearby star-forming cloud. This distance is used to 
estimate the complex mass but, it is important to remark 
that evidence of reddening suggests that d ust material i s 
present up to a distance of about 170 pc ( Franco 1990| ). 
Masses could therefore be underestimated by a factor ^ 3 



1996). I have discovered, here, a sixth cloud which hap- 
pens to be as massive as the Lupus I cloud, ~ 10 4 Mq 
(assuming it is also located at 100 pc). The measured ex- 
tinction for this cloud reaches 4.8 magnitudes. 

The comparison between the 13 CO and the extinction 
map is striking, especially for the Lupus I cloud for which 
each core detected in the molecular observations has a 
counterpart in extinction. Mass estimations can be com- 
pared on condition that the same field is used for both 
maps. Moreover, 13 CO observations are less sensitive than 
B star counts for low extinction. The lower contour in the 
13 CO map corresponds to a visual extinction of ~ 2 magni- 
tudes. Using this iso-extinction contour to define the edge 
of the cloud and the same distance (150 pc) as Tachihara 
et al. ( 1996 ), I find a mass of ~ 1300Mq in agre ement 
with their estimation of 1200M Q . Murphy et al. fll986| ) 
estimate the mass of the whole complex to be ~ 3 10 4 Mq 
using 12 CO observations and a distance of 130 pc. Using 
the same distance, I would obtain 4 1O 4 M (2.3 1O 4 M 
for a 100 pc distance). 



4-8. p Ophiuchus 

Because of the important star formation activity of the 
inner part of the cloud, the IRAS flux at 100^m shows dif- 
ferent structures of those seen in the extinction map. On 
the other hand, the 3 large filaments are present in both 
maps. Even if these regions are complex because several 
stars heat them, the comparison between the far-infrared 
emission and the extinction should allow to derive the 
dust temperature and a 3-dimensional representation of 
the cloud and of the stars involved in the heating. Unfor- 
tunately, uncertainties about the distances for these stars 
are too large (about 15%) and corresponds roughly to the 
cloud size. 

4-9. Orion 



Maddalena et al (1986) have published a large scale CO 
map of Orion and Monoceros R2. Masses derived from CO 
emission are consistent with those I obtain : 1.910 5 Mq 
and 0.86 10 5 Mq for Orion and Monoceros R2, respec- 
tively, from CO data and 3 1O 5 M and 1.21O 5 M from 
the extinction maps. The Orion B maps looks very alike. 
The Orion A maps show a significant difference near the 
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Trapezium (a = 5 h 35 m , 5 = — 5°23') where the young 
stars pollute the star counts. The correlation between CO 
and extinction maps for Monoceros R2 is less striking, 
because of the star forming activity. It is clear that star 
clusters involve an underestimation of the extinction, but 
I would like to stress the point that the heating by a star 
cluster may also destroy the CO molecules. Estimation of 
the column density from star counts and from CO obser- 
vations may, therefore, be substantially underestimated in 
regions such as the Trapezium. 

4-10. Taurus 



Onishi ct al. (1996) have studied the cores in the Taurus 
cloud using a C 18 survey. All of the 40 cores identified 
in their survey are also detected in the extinction map. 



Moreover, Abergel et al. (1994) have shown a strong cor- 
relation between the far-infrared and the 13 CO emission 
in that region. Despite the complexity of the Taurus struc- 
ture (filaments, cores), it is a region, like the Chamaeleon 
complex, located at high galactic latitude (b ~ —16°), 
without complex stellar radiation field. This situation is 
highly favourable to a large scale comparison of CO, far- 
infrared and extinction maps. 

5. Conclusion 

Star counts technique is used since the beginning of the 
century and is still a very powerful way to investigate the 
distribution of solid matter in molecular clouds. Now, with 
the development of digital data, this technique become 
easy to use and can probe much larger areas. For all re- 
gions, we have assumed that all stars were background 
stars. The error resulting from this hypothesis can easily 
be estimated. Equation (|lj) can be written : 

A x = i log (J^J + cste{D ref ) 

where S is the surface which contains nb stars. If 50% 
of the stars are foreground stars, the difference between 
the real extinction and the extinction which assumes all 
background stars is : 

AA X = ilog2 
a 

That corresponds to ~ 0.6 or ~ 1.1 magnitudes of vi- 
sual extinction whether star counts are done using B or R 
band, respectively. Fortunately, most of the clouds are lo- 
cated at small or intermediate distances to the Sun (except 
the Carina at 2500 pc) and this effect is probably small, 
at least for low extinction. The good agreement between 
mass derived from extinction or from CO data argues in 
that favour. 

Stars physically associated to the clouds are more prob- 
lematic because they are located precisely close to the ex- 
tinction cores. Young objects are generally faint in the 



optical band so this problem may be neglected for opti- 
cal counts. This is no longer true with near-infrared data 
for which young objects must be removed before the star 
counts. We obtain extinction map with a spatial resolution 
always adapted to the local densities which are typically 
about ~ 1' for the outer part of cloud and ~ 10' for the 
most extinguished regions where the stellar density be- 
comes very low, i.e. Ss 1000 stars. deg -2 . These maps allow 
the estimation of the total mass of the cloud by extrap- 
olation of the distribution of matter with the extinction. 
It appears that mass concentrated in the regions of low 
extinction represents an important part of the total mass 
of a cloud (1/2 is contained in regions of extinction lower 
than 1 magnitude). 

Extrapolation of the distribution of matter in highly 
extinguished areas is more risky. Star counts method give 
a relation for which masses are underestimated in the 
cores of the cloud for a well understood reason: estima- 
tion of the local density requires to pick up enough stars 
and thus, to use larger area because of the low number 
density. Moreover a characteristic size scale in the distri- 
bution of matter ( Blitz and Williams 1997| ) indicates the 
presence of the lower limit of the fractal cloud structure. 
For this size scale the linear extrapolation used in Fig. [l(] 
also underestimates the real mass and extinction. Despite 
this difficulty, the extrapolated extinction is useful to es- 
timate the saturation level in the extinction map, but it 
is important to keep in mind that extinction can be much 
larger in small cores. 

Finally, the examination of the relation between mass 
and extinction is useful to check what we are measuring. 
The Carina cloud show an aberrant slope (Table |l|) which 
is a strong indication that the map cannot be directly 
interpreted as an extinction map. In that case, the elim- 
ination of reflection nebulae is probably a solution and 
therefore, near-infrared data are requested to investigate 
the extinction. For the Vela and the Serpens cloud, the 
absence of linear part is not understood and reflection is 
not the solution in these regions. 
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